Kinky Behavior in Josephson Junctions 
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We analyze nonperturbatively the behavior of a Josephson junction in which two BCS supercon- 
ductors are coupled through an Anderson impurity. We recover earlier perturbative results which 
found that a 8 = n phase difference is preferred when the impurity is singly occupied and the on-site 
Coulomb interaction is large. We find a novel intermediate phase in which one of 8 — and 8 = n 
is stable while the other is metastable, with the energy E(S) having a kink somewhere in between. 
As a consequence of the kink, the I — V characteristics of the junction are modified at low voltages. 

PACS numbers: 74.50. +r, 73.40.Gk 



I. INTRODUCTION 



Two identical superconductors separated by a bar- 
rier will generate a Josephson current J (6) = J c sin(<5), 
where 6 is the phase difference between the superconduc- 
tors. J c is given by the Ambegaokar-Baratoff formula 
J c = nA/2eR, where A is the superconducting gap and R 
is the normal state resistance of the barrier. This result is 
derived perturbatively, assuming some (spin-conserving) 
tunneling amplitude t for electrons to move across the 
barrier region; the conductance G — 1 /R is then propor- 
tional to \t\ 2 when t is small. In the 1960's, Kulik @] 
showed that tunneling processes which do not conserve 
spin act to decrease J c and potentially drive it negative, 
so that J c oc \t\ 2 — \t sf \ 2 , where t sf is the amplitude for 
spin flip tunneling across the barrier. 

This result was generalized to the case of tunneling 
through a dynamical (Kondo or Anderson) impurity by 
Shiba and Soda Q and by Glazman and Matveev ||. 
The essence of negative Josephson coupling was also elu- 
cidated by Spivak and Kivelson 0, who considered an 
interacting barrier region comprised of a single Anderson 
impurity. Assuming the local Hubbard U is large and the 
impurity is singly occupied, they showed that the fourth 
order process 
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a ring with a single such 7r-junction, time reversal is bro- 
ken and there is a trapped flux of ±/ic/4e In recent 
years, several experiments experiments Q have focussed 
on the existence of 7r-junctions in high T c superconduc- 
tors and have been invoked as evidence for the existence 
of a d x2 _ y2 wave order parameter therein. 

In this paper, we are again concerned with simple s- 
wave superconductors coupled via an Anderson impu- 
rity. We rederive earlier results on 7r-junctions using 
nonperturbative techniques, and we find a crossover re- 
gion, which separates 0- and 7r-junction behavior. This 
crossover regime is characterized by a ground state en- 
ergy E{8) which has local minima at both 5 = and 
5 = tt, and a kink at its maximum, occuring at an in- 
termediate value of 5. We derive a phase diagram as a 
function of — e /A, U/A, and T/A, where T is the bare 
impurity level width generated by virtual hopping onto 
the superconductors. We also discuss the I(V) behavior 
of such a junction within the resistively shunted junction 
(RSJ) model. 



II. EFFECTIVE ACTION 



Consider two superconducting planes each of which is 
connected via electron hopping to an Anderson impurity. 
For electrons living in a fully two-dimensional space, the 
field operator ip( r ) may be decomposed into partial waves 
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reverses the order of spins in the Cooper pair, and hence 
leads to a negative Josephson coupling. For an Ander- 
son impurity with site energy e , the condition for single 
occupancy (for infinitesimal hopping) is U > — e > 0. 

When J c < 0, the most favorable configuration for the 
junction is one in which the phase difference is 6 = tt. In 



where {i/j,(k),ipi(k')} = 27r S(k — k') and where spin 
indices have been suppressed. Within each / sector, the 
kinetic energy can be linearized to s ~ Tiv w {k — k F ). Ex- 
tending the lower limit of k from to — oo is innocuous 
provided we are interested in low energy properties which 
only involve states for which \k — k F \ <C k F . Each I sector, 
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then, yields a single chiral fermion branch on the real line 
(or two chiral fcrmions on a half line) : 
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T l = %v v Jdxipj(x) 0^ - k^j ip^x) . (2) 



Since ip(0) ~ \fk^i/j l=o (0), so only the I = channel is 
involved in connecting the two planes. 

Modeling the superconductors with the BCS Hamilto- 
nian, we make use of 
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FIG. 1. Ground state energy for junctions with e = —2, 
T/A = 1, and (a) U = 2.1, (b) U = 2.6, and (c) U = 2.9. 



which is nonlocal in this representation. It may be ren- 
dered local by means of the canonical transformation 
if>^(x) — > ip^(-x) , which converts up spins into left 
movers. We now integrate out the Fermi fields at all 
points other than x — 0; this means that only the I = 
modes contribute. Defining V4-( T ) = W =( j a {x = 0, r) , 
we obtain the effective action for the point x = 0, 
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^i(uJ m ) RiAuJm) 1pj(Vm) 



(4) 



where 




—iuj r . 



(5) 



Here, aj m = 27t(to + i)T is a fermionic Matsubara fre- 
quency, {ip a ,ip !y } are Grassmann fields, and S is the phase 
of the superconducting condensate (A is real). 

The superconductor is thus reduced the zero- 
dimensional action of equation (|B|). The on-site inter- 
actions of the Anderson impurity and its coupling to the 
superconductors is described by the Hamiltonian ptj 



: E 



{i>Ua + 4^) + £ E 4 c a 



+ Uc \ C \ C l C 1 



(6) 



When t = 0, the ground state of the impurity is singly 
occupied provided U > — e > 0. We couple the Ander- 
son impurity to two superconductors, assumed identical 
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FIG. 2. Josephson current (in units of 2eA/h) for junc- 
tions with e = -2, T/A = 1, and (a) U = 2.1, (b) U = 2.6, 
and (c) U = 2.9. 



in every respect except in their phase. We may then 
integrate out the {ip a ,ip ty } fields on each superconduc- 
tor and decouple the interaction term via a Hubbard- 
Stratonovich (HS) transformation. Neglecting temporal 
fluctuations in the HS field, we have 
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FIG. 3. Occupation number n and magnetization 
m = "f*/U for a junction with e = —2, T/A = 1, and (a, a') 
U = 2.1, (b,b') U = 2.6, and (c,c') U = 2.9. 



where M(uj m ) is the matrix 



-iw m (l + a m ) + e Q + 7 -a m Acos(5/2) 
-a m A cos(5/2) -ia; m (l +a m ) — e + 7 



where a m = a(w m ) = T/y/uj^ + A 2 and T = t 2 /k w v w is 
the (bare) impurity level width. (5 = 8\ — 82 is the phase 
difference between the two superconductors, and e = 
£ + hU. We now make an approximation by performing 
the integral in (Q) using the method of steepest descents. 
That is, we compute the free energy 



F = -T^m([w m (l + a, 



-[a m Acos(i ( 5)] 2 )+^ 7 +e , (8) 



and minimize over 7. We remark that this procedure is 
exact in the large- N limit of an Sp(27V) model in which, 
for each f and J, internal degree of freedom there are N 
flavors a = 1, . . . , JV. The on-site interaction term for 
this model is 
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The free energy in (@) is then the free energy per flavor, 
and e = e Q + (U/2N). 
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FIG. 4. The phase diagram for T/A = 1.0. In the 
phase (which includes all e > as well) ,5 = minimizes the 
energy of the junction and 5 — ir is unstable. In the n phase, 
S — n minimizes the junction energy and S — is unstable. 
In the 0' phase, S = is the stable minimum while S — tt is 
metastable, and vice versa for the tt' phase. 




FIG. 5. The phase diagram for U/A = 12. The phases 0, 
0', n, and 7r' are as in FIG. 0. 



III. SOLUTION OF THE MODEL 

At T = 0, the self-consistent equation for 7 is 

00 

dE 1 fduj 2uj 2 [1 + a(oj)} 2 A(u) 



<9 7 2 



2U 



tt A 2 {Lo) + A-f 2 uo 2 [l + a{uj)f 
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A(to) =lu 2 [1 + a{uj)Y - 7 2 + e 2 + a 2 {oj)A 2 cos 2 (-S) (10) 

Once the solution 7 = 7* is obtained, we compute the 
current J = dE/dS and the impurity occupancy n = 
dE/de , 



J = 



n = 1 



dw A(uj)ol 2 {uj) A 2 sin(<5) 
~ A 2 {uj)+A-f 2 u 2 [l + a{uj)} 2 

00 

duj e A(uj) 

~ A 2 (uj) + 4<y 2 uj 2 [l + a{uj)} 2 



(11) 
(12) 



In the event that there is more than one solution for 7*, 
we choose the solution with the lowest energy. 

In FIGs. [l], H and || we show the mean field solu- 
tion as we tune through a critical region of interaction 
strength U. For low U the junction is unremarkable, 
with J(S) « J c s'm(5) and J c > 0. (The deviation from 
a pure sin(i5) behavior occurs because T/A = 1 is not so 
small.) The energy E(5) exhibits a minimum at 6 = 0, 
and the curve n(S) shows that the Anderson impurity is 
compressible in that the occupancy responds to changes 
in 5. At high U, the behavior is inverted: J c < and 
E(S) has a minimum at S = tt. The occupancy of the 
Anderson impurity is effectively pinned - the impurity is 
incompressible. This is the regime discussed by Spivak 
and Kivelson - for e Q < and U sufficiently large, the im- 
purity maintains single occupancy and the fourth order 
process (in t) which transfers a Cooper pair from one su- 
perconductor to the other through the impurity reverses 
the order of the up and down spins, thereby leading to a 
negative Josephson coupling. 

The intermediate U regime shown in the (b) panels 
of FIGS. |l|, ^, and |3| exhibits unusual behavior. Rather 
than the amplitude of J c smoothly going through zero, 
we find that for S E [0, S c ] the impurity is compressible 
and J(S) > 0, but further increase of 5 results in an in- 
compressible impurity and a reversal of the Josephson 
current. For the parameters given, E(S) has a global 
minimum at 5 = and a local minimum at S = tt; the 
relative positions of the minima will switch at a some- 
what higher value of U. 

A related result was derived in Ref. |l0|, where the 
Josephson tunneling through a ferromagnet was studied. 
The Josephson current was obtained as a sum of partial 
contributions due to Andreev bound states in barrier re- 
gion. Certain bound states were found to produce an 
anomalous phase dependence of the energy, similar to 
that found here. 

The phase diagram for T/A = 1 is shown in FIG. ||. 
In the phases marked and tt, the only stable minima 
in E(5) lie at S = and S = tt, respectively. In the 
intermediate regime, both 6 = and 5 = tt are local 



minima; the global minima then label the phases 0' and 
tt'. We have found that the qualitative shape of this 
phase diagram is valid for all T/A we explored. The 
phase diagram in the (— e /A,r/A) plane is shown in 
FIG. H The fact that the 0' phase persists in an ever- 
narrowing region about e = — ^ U is an artifact of the 
large-TV model, which is unable to describe the Kondo 
physics which sets in below temperatures on the order 
of T K - T exp(-7r|e |/2r) @. The solid lines in FIGS. 
^ and ^, all represent first order transitions, where all 
quantities (n, J, m, etc.) are discontinuous. 

There is a symmetry in our model under e — > — i , 
i.e. — e — * U + e , as is evident in FIG. ||. The phase 
boundaries in the lower left corner all start out linearly 
with the same slope, which may be determined via evalu- 
ating the energy shifts of the impurity states |0) and | f) 
doing perturbation theory in the hopping t. This line is 
determined to 0(T) by Et \ = -Ejj-\, which gives 



F 



U 



ln[ U + A + Vur^] (13) 



4tt y/u 2 - A 2 \U + A — VU 2 - A 2 



This gives a slope dT/d(~ e ) which depends on U/A. 
Our large- ./V theory tends to underestimate this slope. 

How might a cusp in E(S) manifest itself in an exper- 
iment on a single junction? Within the RSJ model, we 
have for I > I c that the average voltage V across the 
junction is 
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V = R 



rds 
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(14) 



For a simple model of the 0' or tt' phase, we take 

{{h/Ae)h5 2 ifO<5<S c 
E„ + {h/Ae)h{8 - tt) 2 if S c < 5 < 2tt - 5 C 
(h/4e)I (6- 2tt) 2 if 2tt - 5 C < S < 2tt 

where Iq, I\, 6 c , and E n are related through continuity 
of E(S) at S = 6 C , the location of the cusp. The function 
I(V) is shown in Fig. ^. The textbook result, when 
E = E (l - cos <5), is I = y/I 2 + (V/R) 2 . In ourmodel, 
the cusp results in a much flatter behavior as V — > 0, 
with 



7(y) = 7 c + A7 c cxp(-V /V) 



(15) 



where I c = max(((5 c /7r) • Iq, (1 — S c /tt) ■ Jij, and A and 
Vo are constants. 

It has been brought to our attention that two conven- 
tional Josephson junctions in series will also yield a kink 
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FIG. 6. I — V diagram for the simple analytical 0' junction 
model (see text). While qualitatively similar to the usual 
/ = ij Ic + (V /R) 2 behavior, here I(V) flattens out much 
more rapidly as V — * 0. The inset shows the energy E(S). 
This I — V behavior also holds for 7r'-junctions. 



in E(S) jy]. If the junctions are identical, then at T = 
each is responsible for half the total phase difference, and 
E(5) — 2Eq(1 — | cos(|<5)|) . However, if each junction is 
resistively shunted, then individual V(I) curves simply 
add, and one does not obtain the behavior in (|l5|). 



IV. CONCLUSIONS 



If single electron tunneling between superconductors 
flips the electron's spin, then a Cooper pair is transferred 
with opposite sign. This result, as first shown by Kulik 
[jjj, can lead to a negative Josephson coupling, in which 
case the ground state of the junction is one in which there 
is a S — 7r phase difference between the two superconduc- 
tors. This phenomenon was subsequently reconsidered 
P,|]j^M], but in a richer and physically realizable con- 
text - tunneling through an Anderson impurity. When 
the tunneling amplitude t onto the impurity is small, the 
ground state is a 7r-junction provided that the impurity 
is singly occupied, i.e. U > —e > 0. We have derived 
the (approximate) nonperturbative phase diagram for the 
junction valid at finite t, and in so doing have discovered 
the existence of two additional phases, denoted 0' and 
7r', in which the ground state energy E(S) has a cusp at 
some 5 £ [0, 7r]. This behavior is manifested in the /(V") 
characteristic of the junction, where for small voltage V 
the current approaches the critical current I c extremely 
rapidly. It is also clear that systems with magnetic im- 
purities should be sensitive to an external magnetic field. 



However, tunneling experiments reported in |12(| demon- 
strated the absence of such dependence up to fields of 6 
T. This could be explained by an antiferromagnetic in- 
teraction between the impurity sites, a possibility we are 
now investigating. 

We gratefully acknowledge conversations with J. 
Hirsch, S. Liu, and D. S. Rokhsar. We are particu- 
larly grateful to S. Kivelson and L. Glazman for sev- 
eral conversations, suggestions, and critical reading of 
the manuscript. 



V. APPENDIX A 



In this appendix we discuss some subtleties associated 
with the Hubbard-Stratonovich transformation for inter- 
acting fermions. Suppose j(t) is distributed according 
to 



P[ 7 (r)] =AAexp ^--L |dr 7 2 (r)j 



(16) 



where J\f normalizes the distribution. Then averaging 
with respect to P[7(r)] one finds, for M(t) a commuting 
variable, 



exp 



dT-f(r) M{t) 



U 







exp | — / drM (r) 



Now consider the simplest Hamiltonian imaginable - a 
single spinlcss fcrmion with 7i a = —fic'c. We compute 
the partition function via a coherent state path integral, 
and since if we let M(t) = c(r) c(r), then from M 2 (t) = 
we obtain the result 



Z a = 



T>[c, c) exp - drc(d T + j(t) - fi) c(r) 



det (a r +7(r)-/i)} 

( f \ 

1 + exp \j3fj. - dr j(t) 



1 + e ,/T e U/2T 



(17) 



which is incorrect - the e u ^ 2T factor should not be 
present in the second term. 

If we consider another trivial model, 



Hb = — m(c|c^ + etc,) + C/c|cJc|C T 



(18) 



and let M(t) = c^c^ — c^Cp we again obtain an incorrect 
result, 
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Z h = l + 2e^ T e u ^ 2T + e 2 ^ T 



(19) 



( 1 - e 7j ) = 1 



These problems do not appear in the discrete time ver- 
sion of the fermionic path integral. The discrete action 
is, for spinless fermions, 



N 



^discrete ^ ^ @i M-ij Cj 



i>3 



where N — (3/e is the number of time slices, and, for Tt a , 



the matrix My is 



/ 1 

-l + erfi 1 

-1 + 672 1 



1 - e lN \ 





V 



-l + e7iv-i 1 / 



where 7„ = 7 n — /i. We must now integrate out the Grass- 
mann variables and then average over the distributions 



for each time slice j. But now we find 

N 

det M = 1 + Y[ (1 - Hj + £M) 



(21) 



and indeed 



Z a = (detAf( 7l ,...,7 JV )) {7j 

= l + (l + e M ) iV 
= 1 + e^ T . 



(22) 



whereas if we write 1 — 67 sa exp(— €7) and then average, 



(exp(-£7j)) = ex P(2 eC/ ) 



(20) which gives an incorrect contribution at 0(e)- 



The remedy for this frustrating problem is to perform 
a shift 



(25) 



in the continuum case. This rescues the correct results 
from both equations (|j]) and (jl^) , and may be formally 
obtained by point splitting M 2 (t) -> Af(r+) M{t~). 



VI. APPENDIX B 



To demonstrate the reliability of the formalism and 
approximations used in this paper, we consider the toy 
model 



/i(c|c| + ejejj + A(clcl + c,cS) + Uc^c^c^c 



t„t 



J J, 



Once again we break up the interaction term with a 
Hubbard-Stratonovich transformation, careful to substi- 
tute All the time dependence of the field 7(1-) 
can be gauged away via the transformation 



c(t) -> e~ 5(r) c(t) 



1 



P 



9r9(r) =7(t)-~ ydr' 7 (r') 




(26) 



When we apply the discrete time path integral to Tib, 
we find 

AT JV 

det M T det M x = 1 + JJ (1 - £7, + e/x) + J| (1 + e7j + e/x) 
3=1 j'=i 



Note that <?(/?) = 5(0), so the antiperiodic boundary con- 
ditions c(f3) — — c(0) are not affected by this transforma- 
tion. Thus. 



1 



d7e -^/2UT e (fi+^/T 



N 



Hiil+erf-e^) (23) 
3=1 



v / 2nUT 

— OO 

x det ( iL A 



3 r + 7 + /i 



(27) 



and we find, correctly, 

Z b = l + 2e^ T + e^- l/ )/ T . (24) 



The essence of the problem is that in the discrete case 
we have 



where now 7 is time-independent, denoting simply the 
average of what we have until now been calling the func- 
tion 7(7"). Now if {Ao, • . • , A3} are constants, then 

det (<9 r + A + A • 3) = 1 + 2 e" Ao/T cosh ^ + e ~ 2A »/ T . 



(28) 
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FIG. 7. Occupancy n versus bare chemical potential /x for 
the toy model TL C computed in the steepest descents approx- 
imation (solid) and compared with exact results (dashed). 



Performing the integral over 7, one correctly obtains the 
partition function which may be calculated by more ele- 
mentary means. 

The steepest descents approximation (SDA) to the in- 
tegral over 7 gives 



2 
u 



1 - e' 2 ^ T 



1 + 2 e-T/r cosh ^ 2 T +A2 + e- 2 ^ T 



(29) 



to isolated discontinuous changes of 7 as fx, U, and A are 
varied, which are reflected in the behavior of e.g. 
The exact solution behaves smoothly, but for sufficiently 
low temperatures the agreement is arbitrarily good. 
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which is obtained by extremizing 

,2 



F = 



7 
2U 



7-A-Tln(l + e -W^ X2 ") /TN 
TlJl + e-W^Ai)/^) (30) 



with respect to 7. 7 = is always a solution, and for a 
range of parameters it is possible to have more than one 
solution to (|9|) , in which case we choose the solution with 
the lowest free energy. Once a solution for 7 is obtained, 
we evaluate the particle number 

2e -7/T smh V£pI 



(31) 



y/jl 2 + A 2 1 + 2e - 7 /T 



cosh j/£^! +e -2 7 /T 



Comparisons of the SDA with the exact results are shown 
in FIG. 0. For the Hamiltonian TL C , the SDA is exact in 
the limit of zero temperatures. At finite T, the SDA leads 
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